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the static deformation field produced by point forces and 
point force couples embedded in two Poissonian half- 
spaces in welded contact. Kumari et al. (1992) 
generalized the results of Heaton and Heaton (1989) 
and they obtained the expressions for the displacements 
for arbitrary Poisson's ratio. Furthermore, they obtained 
closed-form analytical expression for the stresses as 
well. Sharma et al. (1991) obtained closed-form 
analytical expressions for the displacement and stresses 
at any point of either of the two homogeneous, isotropic, 
perfectly elastic half-spaces in welded contact due to 
very long strike-slip dislocations. 



I. ABSTRACT 

Closed-form analytic expressions for the 
Airy stress function, stresses and displacements for 
half-spaces in smooth contact by applying suitable 
boundary conditions at the interface are first 
obtained. These expressions are then integrated 
analytically to derive the Airy stress function for a 
vertical tensile fault of arbitrary dip and finite width. 
Closed-form analytical expressions for the 
displacements and stresses follow immediately from 
the Airy stress function. The variation of the 
stresses with the horizontal distance from the fault 
and with the depth studied numerically. The results 
obtained in smooth contact are compared with that 
of welded contact. 

KEYWORDS Deformation, Dyke injection. Vertical 
tensile fault. Smooth contact. Welded contact 



II. INTRODUCTION 

Rongved (1955) obtained closed-form algebraic 
expressions for the Neuber-Papkovich displacement 
potentials for an arbitrary point force acting in an infinite 
medium consisting of two elastic half-spaces in welded 
contact. Dundurs and Hetenyi (1965) obtained these 
functions when the half-spaces are in smooth contact 
and also obtained the corresponding displacement field. 
Moreover, they compared the traction for a smooth 
interface with that for a perfect bond between the two 
semi-infinite solids. Heaton and Heaton (1989) used the 
expressions for the Neuber-Papkovich functions derived 
by Rongved (1955) to obtain analytical expressions for 



Singh and Garg (1986) obtained the integral 
expression for the Airy stress function in an unbounded 
medium due to various two-dimensional seismic 
sources. Beginning with these expressions. Rani et al. 
(1991) obtained the integral expressions and hence 
closed form analytical expressions for the Airy stress 
function, displacements and stresses in a 
homogeneous, isotropic, perfectly elastic half-space due 
to various two-dimensional sources by applying the 
traction-free boundary conditions at the surface of the 
half-space. Singh and Rani (1991) obtained closed-form 
analytical expressions for the displacements and 
stresses at any point of a two-phase medium consisting 
of a homogeneous, isotropic, perfectly elastic half-space 
in welded contact with a homogeneous, orthotropic, 
perfectly elastic half-space caused by two-dimensional 
seismic sources located in the isotropic half-space. Singh 
et al. (1992) followed a similar procedure to obtain the 
closed-form analytical expression for the displacements 
and stresses at any point of either of two homogeneous, 
isotropic, perfectly elastic half-spaces in welded contact 
due to two-dimensional sources. 



17 



www.ijitce.co.uk 



International Journal of Innovative Technology and Creative Engineering (ISSIM:2045-871 1) 

VOL.3 No.1 January 2013 



Singh and Singh (2000) derived closed-form 
analytical expressions for the displacements and 
stresses at an arbitrary point of the half-space caused 
by a long vertical tensile fault of finite width. Singh et al. 
(2002) extended the results of Singh (2000) and 
obtained the closed-form analytical expression for the 
subsurface stresses and displacements caused by a 
long inclined tensile fault buried in a homogeneous 
isotropic half-space. The effect of depth of the upper 
edge of the fault and the dip angle on the deformation 
field has been examined. Kumar et al. (2005) derived 
analytical expressions for the displacements and stress 
fields at any point of the two homogeneous, isotropic 
elastic half-spaces in welded contact caused by a long 
tensile fault of arbitrary dip and finite width. 

In this paper, we obtain closed-form analytical 
expressions for the displacements and the stresses at 
any point of either of the two homogeneous, isotropic, 
perfectly elastic half-spaces in smooth contact caused 
by various two-dimensional sources embedded in one of 
the half-spaces. The results are valid for arbitrary 
Poisson's ratio. Only plane strain case is considered. 
The stresses for the smooth interface case are 
compared with the corresponding stresses for the 
welded interface. 

We begin with the integral expressions for the 
Airy stress function in an unbounded medium given by 
Singh and Garg (1986) to obtain the integral 
expressions for the Airy stress functions stresses and 
displacements for the two half-spaces in smooth contact 
by applying suitable boundary conditions at the 
interface. The integrals are then evaluated analytically, 
obtaining closed-form expressions for the stresses at 
any point of either of the two half-spaces caused by 
vertical tensile fault. For numerical calculations, half 
spaces are assumed Poissionian. 



III. THEORY 

Let the Cartesian coordinates be denoted by 
(1^x2,13) with xg-axis vertically downwards. Consider 
a two-dimensional approximation in which the 
displacement components u^, wg, u-^are independent of 

A-^ so that 3/3t^ = 0. Under this assumption, the plain 

strain problem (» j = 0)can be solved in terms of the Airy 
stress function u such that 



d u 

^22 - ^TX ' ^33 ■ 



d u 



d u 
d.v2av3 



0 
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where r are the components of stress and 



72 = 1! + ^ 

~ 3 2 ^ 2 

dv2 dv3 



We have (Sokolnikoff,1956,Sec.71) 
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where a = 



U+2n) 2(1- (t) 
We consider two homogeneous, isotropic, perfectly 
elastic half-spaces which are in smooth contact along 

the plane 13 = 0. The upper half-space (a3( 0) is called 
medium I and the lower half-space (xg) 0), medium II, 

with elastic constants /Ij, //j and ^,1^-2 respectively 

(Figure 1). The superscript (1) denotes the quantities 
related to medium I and the superscript (2) denotes 
those related to medium II. 

Let there be line source parallel to xj-axis. As 
shown by Singh and Garg (1986), the Airy stress 
function u q\q\ a line source parallel to the x^axis 
passing through the point (0,0, a) in an infinite medium 
can be expressed in the form 

-A I) sin/f.i 



U , = \ |j(in-l-M ,Vf|x 



(4) 



r'i 



+(^o+eo4^"3-''|) 

where the source coefficients z,q, m q, /'q and Qgare 

independent of the variable * .For a line source parallel 

to the X J axis passing through the point (0,0, a) of 

medium II, the expression for the Airy stress functions in 
the two half-spaces are of the form 

u = j[(ij-i-A/ jA-x3)sini-A2 

0 

+(^l+Oi/fA3)cOSAA2j A-"\-*'3ciA- , 



(5) 
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(6) 



0 

p 2 + e 2* -^^ 3 ) CO s/f A- 2] i- " 3 d- , 
where u ^ is given in equation (4) and tlie 
unknowns!^, M ^, and i^, m ^, p^, g ^are to be 

determined from the boundary conditions. 
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Fig. 1 . Two half-spaces in smootli contact with a line source in the 
lower half-space at (0, 0, h) 

Since the interface is vg=0, then smooth contact implies 

that Uganda are continuous but vanishes at the 

interface. Since the half-spaces are assumed to be in 
smooth contact along the plane v,^=0, the boundary 

conditions at the interface v„=0 are 
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Due to discontinuity in horizontal displacement « 2 

at v,^ = 0, there may be a slip between the two media. 

From equations (1) - (6), we get the following integral 
expressions for the stresses and the displacements for 
the two media: 



M I 
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(2) 
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It is noticed from appendix that the 
coefficients Iq, M Q, andQphave different values 

forx3<A; let i~,M~,p~and g be the values of 

Lq, M Q, Pq and Qq respectively, valid forA3<A. Then r*g'=-c^ 

from equations (9) to (17) and the boundary conditions 
(7), we obtain 
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where 



(-- ;^A),i^^ = .v-2+(.x3-0^ S^=-2+(-3 + 0^ 

The expressions given in equations (18)-(29) 
are in terms of the source coefficients L„, P„ and Q„. 
The results for the Airy stress functions and the stresses 
in the two half-spaces for a vertical tensile source on 
putting the values of the source coefficients from given 
appendix are as follows: 
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IV. NUMERICAL RESULTS 



We want to compare the stresses for the smooth 
interface case with the corresponding stresses for the 
welded interface. For the numerical computations, we 
assume that the half-spaces are Poissonian(A =//Jso 

2 



that ctTi = cir2 ■ 



3 



^ l+2p 



-(ci+fli) 



We define the following dimensionless 
quantities!' = — , Z = — , where h is the distance of the 

A A 

line source from the interface. Let the dimensionless 
stresses for the welded interface(Kumar et.al 2005) be 
denoted by p **^^ and the dimensionless stresses 



for the smooth interface be denoted by p^^\pi^^ . Then 



(1) ^(1) •(!) •(!) 

P = T , P = T 

'> . . 'J , , '-' 



(38) 



From equations (31)-(33) and (35)-(37),we get the 
following expressions for the dimensionless stresses for 
the two cases for a vertical tensile line source: 
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Following are the expressions for welded contact 
(Kumar et.al 2005) for the dimensionless stresses for 
the two cases for a vertical tensile line source: 
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16c oZ (z + Dd- ^ + 1+ Z ) 8c oF ^(z ^ - 1+ 2z ) 



96c2zr^(z + lf 



(50) 



V. DISCUSSION 



Figures 2(a)-2(d), 3(a)-3(d) and 4(a)-4(d) show the 
variation of the dimensionless stresses P^^, P23and P^^ai 
the interface with the horizontal distance from the fault 
and the figures 5(a)-5(d), 6(a)-6(d) and 7(a)-7(d) show 
the variation of the dimensionless stresses P^^, P^^ and 
P33 at the interface with the depth from the fault for the 



rigidity contrast 



M2 



1 3 

-, 1 -, 5, respectively for the two 

5 2 



cases. In all the figures (from 2(a)-7(d)),f^)^^and Pif^ are 
dimensionless stresses in the case when contact of half- 
space is smooth, whereas /^J'^'and /^j''" denote 

stresses when contact of half-space is welded. The 
superscript denotes the medium. Figure 2(a) shows 
variation of normal stress P^^ for smooth contact and 
in case of welded contact with distance from the fault for 



M2 5 



In the succeeding figures i.e. 2(b) to 2(d), we 



increase the rigidity ratio to 1,—, 5 respectively. It is 

noticed that near the origin, the values of stress 
components pS^- and pSP lies between P*}^- andP-t'^-. 
As we move away from the fault, this difference 
becomes less and approaches to zero for infinitely large 
distances from the fault. Figures 3(a) to 3(d) show 
variation of stress /^gand ^2*3 with distance from the fault 

for the rigidity contrast — = -, 1, — and Srespectively. 

//2 5 2 

As the value of rigidity contrast increases the value of 
stress components P^}^- and F^}^- decreases in welded 
contact. It is also observed in figures from 2(a) to 4(d) 
that initially, the value of stress is larger in medium I as 
compared to medium II for smooth contact, components 
FjHs, P^and P^. 

Variation of stresses with depth from the fault is 
shown in figures 5(a) to 7(d). Here, too, variation is 
examined for same rigidity ratios as explained above. It 



is observed that initially the values of stresses P-~ar\6 
Pg3 are nearly zero. 

VI. Conclusion 

In both cases i.e in welded contact and 
smooth contact, stresses tend to zero as x^ and x, 
approaches to infinity respectively. We notice that the 
magnitude of stresses for smooth contact model are 
greater than the magnitude of stresses for welded 
contact model up to some epicentral distances and then 
decreases rapidly. Also, in case of smooth contact, the 
behavior of stresses at the interface for the two half- 
spaces is opposite to each other i.e., when it is 
increasing in medium I, it is decreasing for medium II 
and vice versa. We observe that a change in value of 



alters the magnitude of the stresses. 



p(i) 
^22 

p(2) 
'^22 

p"(l) 
22 

p"(2) 
^22 
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Fig. 2(a) Variation of tine dimensionless normal stress P^^ with 
distance from tfie fault for jUi/jU2 = y5 for smootfi contact and 
welded contact. 
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Fig. 2(b) Variation of tfie dimensionless normal stress P^^ witfi distance 
from tfie fault for ,Uj///2 = l for smootfi contact and 
welded contact. 
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Fig. 2(c) Variation of tine dimensionless normal stress P,, with Fig. 3(b) Variation of the dimensionless shear stress P,, with 
distance from the fault for = 3/2 for smooth contact and distance from the fault for //1///2 = 1 for smooth contact and 

welded contact. welded contact. 
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Fig. 2(d) Variation of the dimensionless normal stress P^^ with distance 
from the fault for //j/ jU2 = 5 for smooth contact and welded contact. 
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Fig. 3(a) Variation of the dimensionless shear stress P^, with distance 
from the fault for //j///2 = l/5 for smooth contact and welded 
contact. 
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Fig. 3(c) Variation of the dimensionless shear stress P^j with distance 
from the fault for jUi/ fJ^ = 3/2 for smooth contact and welded 
contact. 
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Fig. 3(d) Variation of the dimensionless shear stress P^j with 
distance from the fault for //1///2 = 5 for smooth contact and 
welded contact. 
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_ p(2) 



Fig. 4(a) Variation of tine dimensionless normal stress P33 witin distance 
from tfie fault for H^j //2 = 1/5 for smootfi contact and welded 
contact. 
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Fig. 4(d) Variation of the dimensionless normal stress P33 with 
distance from the fault for fA-^j //2 = 5 for smooth contact and 
welded contact. 
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Fig. 4(b Variation of the dimensionless normal stress P33 with distance 
from the fault for //^///2 = 1 for smooth contact and welded contact. 
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Fig. 4(c) Variation of the dimensionless normal stress P33 with 
distance from the fault for jA-^j IJ.2 = 'il2 for smooth contact 
and welded contact. 
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Fig. 5(a) Variation of the dimensionless normal stress P^^ with 
depth from the fault for JU-yl /U2 = ^/ ^ for smooth contact and 
welded contact. 
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Fig. 5(b) Variation of the dimensionless normal stress P^^ with 
depth from the fault for JUij Hi — ^ for smooth contact and 
welded contact. 
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Fig. 5(c) Variation of tine dimensionless normal stress P^^ witli 
deptli from tfie fault for //j///2 = 3/2 for smootfi contact 
and welded contact. 
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Fig. 5(d) Variation of tfie dimensionless normal stress P^^ witfi 
deptfi from tfie fault for /Ji/ //^ = 5 for smootfi contact and 
welded contact. 
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Fig. 6(a) Variation of tfie dimensionless sfiear stress P^j witfi 
deptfi from tfie fault for //j///2 = 1/5 for smooth contact 
and welded contact 
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Fig. 6(b) Variation of the dimensionless shear stress P^j with 
depth from the fault for ^j///^ = 1 for smooth contact and 
welded contact. 
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Fig. 6(c) Variation of the dimensionless shear stress P^j with 
depth from the fault for f.li/ //^ = 3/2 for smooth contact 
and welded contact. 
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Fig. 6(d) Variation of the dimensionless shear stress P^j with 
depth from the fault for /Ui/ H'l — ^ for smooth contact and 
welded contact. 
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Fig. 7(a) Variation of tine dimensionless normal stress P33 with 
deptli from tfie fault for //j/ //^ = 1/5 for smootfi contact 
and welded contact. 
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Fig. 7(d) Variation of tfie dimensionless normal stress P^, with 
depth from the fault for /^i///2 = 5 for smooth contact and 
welded contact. 



APPENDIX 
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Vertical tensile fault 
Horizontal tensile fault 
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Fig. 7(b) Variation of the dimensionless normal stress P33 with 
depth from the fault for //j///2 = l for smooth contact and 
welded contact. 
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The upper sing is for .v.^ > h .the lower sign is for <h,b is the 

magnitude of the displacement dislocation and d is the width of the 
line fault. 
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Fig. 7(c) Variation of the dimensionless normal stress P33 with 
depth from the fault for fl^j //^ = 3/2 for smooth contact 
and welded contact. 
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